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The parameters of the Galactic spiral wave are re-determined using a modified periodogram (spectral) analysis of the 
galactocentric radial velocities of 58 masers with known trigonometric parallaxes, proper motions, and line-of-site veloc- 
ities. The masers span a wide range of galactocentric distances, 3 <R< 14 kpc, which, combined with a large scatter of 
position angles 6 of these objects in the Galactic plane XY, required an accurate account of logarithmic dependence of 
spiral-wave perturbations on both galactocentric distance and position angle. A periodic signal was detected corresponding 
to the spiral density wave with the wavelength A = 2.4 ± 0.4 kpc, peak velocity of wave perturbations /r = 7.5 ±1.5 
km s~\ the phase of the Sun in the density wave x© = —160 ± 15°, and the pitch angle of —5.5 ± 1°. 
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1 Introduction 

Spectral analysis of residual velocities of different young 
galactic objects (HI clouds, OB stars, open star cluster 
younger than 50 Myr, masers) tracing the Galaxy spiral 
arms has been fulfilled, for example, by Clemens (1985), 
Bobylev, Bajkova & Stepanishchev (2008), Bobylev & Ba- 
jkova (2010). As a result there were determined the fol- 
lowing parameters of the spiral density wave subject to the 
theory by Lin & Shu (1964): amplitude and wavelength 
of the perturbations, evoked by the spiral wave, pitch an- 
gle, phase of the Sun in the spiral wave. Nowadays galac- 
tic masers having high-precision trigonometric parallaxes, 
line-of-sight velocities and proper motions (Reid et al. 2009, 
Rygl et al. 2010) are of great interest. 

The previous spectral analysis represents the simplest 
periodogram analysis of velocity perturbations based on 
conventional Fourier transform, what can be considered 
only as the first approximation of exact spectral analysis and 
can be applied adequately only in the case of small range of 
galactocentric distances (2-3 kpc). But analysis of modern 
data on galactic masers which are located in wide range of 
galactocentric distances (3 <R< 14 kpc) requires elabo- 
ration of more correct tools of spectral analysis accounting 
both logarithmic dependence from galactocentric distances 
and position angles of objects. For a detailed description of 
the new method of spectral analysis of velocity residuals see 
Bajkova & Bobylev (2012). 



Our first study (Bobylev & Bajkova 2010) was based 
on an analysis of radial galactocentric velocities of only 28 
Galactic masers. The second one (Bajkova & Bobylev 2012) 
dealt with 44 masers. Currently, high-precision VLBI mea- 
surements of parallaxes, line-of-site velocities, and proper 
motions are available for 58 Galactic masers, which is of 
great interest for our task. The aim of this present study is to 
re-determine the spiral density wave parameters by applying 
the recently proposed algorithm (Bajkova $ Bobylev 2012) 
and new ones, described below, to more extensive data se- 



2 Basic relations 

The velocity perturbations of Galactic objects produced by 
a spiral density wave (Lin & Shu 1964) are described by the 
relations 



Vr = -f R cosx, 
AV e = fesinx, 



(1) 

(2) 



where 
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X = m[cot(») HR/R ) -6}+ XQ (3) 

is the phase of the spiral density wave; m is the number of 
spiral arms; i is the pitch angle; x© is the phase of the Sun 
in the spiral density wave (Rohlfs 1977); R is the galacto- 
centric distance of the Sun; 9 is the object's position angle: 
tan 8 — y/(R — x), where x, y are the Galactic heliocen- 
tric rectangular coordinates of the object; Jr and fg are the 
amplitudes of the radial and tangential perturbation compo- 
nents, respectively; R is the distance of the object from the 
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Galactic rotation axis, which is calculated using the helio- 
centric distance r = 1/ir: 



where 



R z = r cos b - 2R r cos b cos I + Ri 



(4) 



where I and b are the Galactic longitude and latitude of the 
object, respectively. 

Equation (01 for the phase can be expressed in terms of 
the perturbation wavelength A, which is equal to the dis- 
tance between the neighboring spiral arms along the Galac- 
tic radius vector. The following relation is valid: 

— - — = mcot(i). (5) 

A 
Equation (0 will then take the form 

2irR 



X = 



X 



■\n(R/R )-m8 + XQ- 



(6) 



The question of determining the residual velocities is 
considered below. To goal of our spectral analysis of the 
series of measured velocities V Rn , AVg n n = 1,2, ... ,N, 
where N is the number of objects, is to extract the peri- 
odicity in accordance with model (Q~|l-(0 describing a spi- 
ral density wave with parameters f R ,f$, X, and \o- If tne 
wavelength A is known, then the pitch angle i is easy to de- 
termine from Eq. (0 by specifying the number of arms m. 
Here, we adopt a two-armed model, i.e., m = 2. 

3 Methods 

3.1 Fourier transform-based analysis 

Let us represent series of velocity perturbations of galactic 
objects evoked by a spiral density wave (Q]i-(|2]l in the most 
general, complex form: 

V n = V Rn +jAV 6n , (7) 

where j = y/—l, n is a number of an object (n — 1, ..., N). 

A periodogram analysis, which we consider here, re- 
quires calculation of power spectrum of series (0 expanded 
over orthogonal harmonic functions 

exp[-j ln(R n /R ) + jmd n ] 

Afc 

in accordance with expression (0 for the phase. 

A complex spectrum of our series is: 



Vx k 



1 N 

71 = 1 



V x f + jHr = - > ; (V Rn + jAV 9n ) x (8) 



x exp[-j— — ^ bi(Rn/R ) + jmd n ], 
Xk 
where the upper indices Re and Im designate real and 

imaginary spectrum parts respectively. 

Let us reduce the latter expression to a standard discrete 

Fourier transform in the following way: 



1 N 
^ = n zJ( Fr ™ +3& V eJ x exp(jm0 n ) 

x exp[-j— — ^ ln(R n /R )] ■■ 



1 N 9ttR 

-^v^M-j^HRn/Ro)], 



(9) 



rRe 



(10) 



V = V" e + v lm = 

= [Vji n cos(m0 n ) — AVg n s'm(m6 n )] + 
+j[Vn n sm(m9 n ) + AVg n cos(m9 n )]. 

And, finally, making the following change of variables 

R n = R \n{R n /R ), (11) 

we obtain standard Fourier transform of a new series V n 
( [Tol l, determined in point set R n : 



N i 

f> ! \"^ t/ r • 2nR r , 

n—1 



(12) 



The periodogram |VaJ 2 is subject to further analysis. 
The peak of the periodogram determines the sought-for pe- 
riodicity. The coordinate of the peak gives the wavelength A 
and, respectively, pith angle i (see Eq.©). Relation between 
a peak value of the periodogram S pea k an d perturbation am- 
plitudes fji and f$ is expressed as follows: 

f 2 R + fe 2 = 2xS peak . (13) 

It is necessary to note that the spectral analysis of com- 
plex series (f7]l allows to determine A (or pitch angle i) and 
phase of the Sun x©, but does not allow to estimate the am- 
plitudes fu and fg separately. To determine them it is nec- 
essary to analyze radial {Vr„} and tangential {AVe„} ve- 
locity perturbations independently, as it has been shown by 
Bajkova & Bobylev (2012). Here we show how amplitudes 
of perturbations can be found if A and \q are known (for 
example, from previous complex analysis). We consider se- 
ries {Vr„} (by analogy the same algorithm can be applied 
to series {AV$ n }). 

Let us represent Eq. (0 as 

m6, (14) 



X = Xi 
where 

27ri?. 



Xi 



X 



]n(R/R ) + XQ . 



(15) 



Substituting (TT4b into Eq. (Q3 for the perturbations at the 
nth point and performing standard trigonometric transfor- 
mations, we will obtain 

V Rn = -/rcos(xi„ - m6 n ) = (16) 

= -IrCOSXi„ cos m6 n - /^sinxu sinm0„ = 

= — /r cos %i ra (cos m^ n + tanxi n smm6 n ). 

Let us designate 

Vr = -/ r cosxi, (17) 

Owing to the substitution (TFTT i. it then follows from ( TTol l that 

Vr„ = V Rn (cos m9 n + tan xi n sin m9 n ). (18) 

Substituting known values of A and x© we can form from 
Eq. ( fT8l a new data series 

V Rn = V Rn /(cosm9 n + tanxi„ sinm9 n ). (19) 
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Again, using the substitution (fTTT i, we obtain a standard 
Fourier transform: 



n=l 



A, 



(20) 



In this case relation between a peak value of the peri- 
odogram S pea k and perturbation amplitudes Jr is as fol- 
lows: 



JR — 4 X bpeak- 



(21) 



Note, that a separate periodogram analysis based on op- 
erations ([T4li-(l20i> can be realized as an iterative process of 
seeking for unknowns A, x©> and Jr under optimization of 
some specific signal extraction quality criterium (Bajkova 
& Bobylev, 2012). 

For numerical realization of Fourier transform (ITZl or 
( f20l > using fast Fourier transform (FFT) algorithms it is nec- 
essary to determine data V n (R )(n — 1, . . . , N) on discrete 
grid I — 1, . . . , K — 2 a , where a is integer, positive, N < 
K; An is a discrete space. Coordinates of data are deter- 
mined as follows: l n = [(R n + \mm{R k }\ k= i i ,,, iN )/A R ] + 
l,n — 1, ..., N, where [a] denotes an integer part of a. The 
sequence determined is considered as a periodical one with 
the period D = K X Ar, Obviously, the values of the K— 
point sequence are taken to be zero in the pixels into which 
no data fall. 



In our case, the reconstruction problem assumes find- 
ing the minimum of the following generalized entropy func- 
tional: 



E=J2V k 

k=Ki 



Re+ HaV k Re+ ) + V k Re - \n(aV Re - 



(24) 



-V, 



Im+ \n(aV k Irn+ ) + V k Im ~ HaV k Im - 



E 

71=1 



Im\2 



{riTY + K™) 



2al 



where the sought-for variables V k Re and V k Im are repre- 
sented as the difference of the positive and negative parts: 

V k Re = V k Re+ - V k Re ~ and V k m = v k m+ - Y k m ~ re- 
spectively; in this case, V Re+ , V Re ~ , V k Im+ , V k Im ~ > 0, 
a > 1 is the real-valued parameter responsible for the 
separation of the positive and negative parts of the sought- 
for variables with the required accuracy (in our case, we 
adopted a = 10 6 ), K\ and Ki are the a priori known lower 
and upper localization boundaries of the sought-for finite 
spectrum, r/ Re and 77^™ are a real and an imaginary parts re- 
spectively of the measurement error of the nth value of the 
series V n that obey a random law with a normal distribution 
with a zero mean and dispersion a n . 

The constraints ( 1221 on the unknowns, with accounting 
measurement errors, can be rewritten as 



3.2 The GMEM-based analysis 

So far we have considered the simplest method of peri- 
odogram analysis based on linear Fourier transform. In the 
case where the data series are irregular, i.e., there are large 
gaps, the signal spectrum is distorted by large side lobes 
and it becomes difficult to distinguish the spectral compo- 
nent of the signal from spurious peaks. In this case, it may 
turn out to be useful to apply nonlinear methods of spec- 
trum reconstruction from the available data. This problem 
is fundamentally resolvable if the sought for signal has a 
finite spectrum. Since our problem belongs to the class of 
problems on the extraction of polyharmonic functions from 
noise, we assume that this condition is met. 

Here we propose a complex spectrum reconstruction al- 
gorithm based on well-known maximum entropy method 
(MEM). Since the spectrum is described by a complex- 
valued function, we apply a generalized form of MEM 
(GMEM) proposed and described in detail by Bajkova 
(1992) and Frieden & Bajkova (1994). 

The spectrum Vk — V Re + jV k Im and the data V n are 
related by the inverse Fourier transform: 

27r(*-l)(l„-l), 



k=K2 



K 

E 



(vr+jvn^u- 



K 



= V n . (22) 



Note that wavelength \ k and spatial frequency (k — 1) 
are related as follows: 
D 



Afe — 



fc-1 



(23) 



£(tt 



Re+ 



V, 



Re--. 



■W' 



v, 



Im — 



k=Kl 



x exp(j 



2ir(k-l)(l T 



>l 



Re 



m 



Im 



)) x (25) 



K. 



We can see from (l24l i. that the functional to be mini- 
mized consists of five terms, the first four ones are total en- 
tropy of the sought-for solution, the last one is \ 2 measure 
of deviation between data and solution. Thus, the GMEM 
algorithm seeks for solutions not only for the spectrum un- 
knowns, but also for measurement errors (rj Re + jrj^ m ). 
Therefore we can expect effective suppression of noise 
caused not only by non-uniformity of series but also by 
measurement errors. The optimization of functional (l24l un- 
der conditions ( |25] | can be done numerically using any gra- 
dient method. We used a steepest-descent method. 

4 Data 

Bobylev & Bajkova (2012) analyzed a sample of 44 masers 
with known high-precision trigonometric parallaxes, line- 
of-site velocities and proper motions. Since then the num- 
ber of masers with such measurements has increased con- 
siderably. Now data on 58 masers are available in literature. 
Table 1 lists the input data on 58 masers associated with the 
youngest Galactic stellar objects (protostar objects of dif- 
ferent mass, very massive supergiants, or T Tau stars). The 
references to majority of original data can be found in Ba- 
jkova & Bobylev (2012). 
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Table 1 Data on masers 



2 4 6 
Y, kDC 



8 10 



Fig. 1 Coordinates of masers in the XY Galactic plane 
(the Sun is located at the center of coordinate system). 

The input data, namely, trigonometric parallaxes and 
proper motions, were obtained by several groups using long 
time radio-interferometric observations carried out within 
the framework of different projects. One of them — the 
Japanese project VERA (VLBI Exploration of Radio As- 
trometry) (Honma et al. 2007) — is dedicated to observa- 
tion of H2O and SiO masers at 22 and 43 GHz, respectively. 
Note that higher observing frequency results in higher reso- 
lution and more accurate data. Methanol (CH3OH) masers 
were observed at 12 GHz (VLB A, NRAO) and 8.4-GHz 
continuum radio-interferometric observations of radio stars 
were carried out with the same aim (Reid et al. 2009). 

The line-of-site velocities V r (LSR) listed in Table 1 
were determined with respect to the Local Standard of Rest 
by different authors from radio observations in CO emission 
lines. The parallaxes were determined, on average, with a 
relative error of u-^/tt ~ 5%, and only in three regions the 
error exceeds the mean level. These are IRAS 16293-2422 
(ct^/tt = 19%), G 23.43-0.20 (ct^/tt = 18%) and W 48 
(a^/n = 14%). 

Fig. Q] shows the space distribution of masers projected 
onto the Galactic XY plane. The objects can be seen to be 
widely scattered along the x, y coordinates, implying a large 
scatter of position angles. Hence the spectrum analysis algo- 
rithm has to be constructed in order to correctly extract pe- 
riodic signal from velocity perturbations traced by masers. 

5 Results 

First, we re-determined the parameters of the Galactic ro- 
tation curve using the data for 58 masers. The method 
employed is based on the well-known Bottlinger formu- 
lae (Ogorodnikov, 1965), where the angular velocity of 
Galactic rotation is expanded into a series up to 2-nd or- 
der terms in t/Rq (Bobylev & Bajkova 2010). We adopt 
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V Rn = -U n cos 9 n + (Vq + 14) sin 6 n , 



(27) 



4 6 8 10 12 

R , kpc 

Fig. 2 Masers galactocentric radial velocities vs galacto- 
centric distances R. 



Ra = 8 kpc and infer the following components of the pe- 
culiar solar velocity: (Uq,V q ,W q ) = (7.4, 16.6, 8.53) ± 
(1.0, 0.8, 0.5) km c _1 , and the following parameters of the 
Galactic rotation curve: flo = — 29.3 ±0.6, km c _1 kpc -1 , 
Q' = +4.2 ± 0.1 km c- 1 kpc ~ 2 , ft' ' = -0.85 ± 0.03, km 
c _1 kpc ~ 3 . The linear Galactic rotation velocity at R = R 
then is equal to: V = \RqC1 \ = 234 ± 5 km c _1 . 

There is good agreement of our results with the re- 
sults of analyzing masers by different authors. Based 
on a sample of 18 masers, McMillan & Binney (2010) 
showed that Q,q lying within the range 29.9 — 31.6 km 
s _1 kpc -1 at various Rq was determined most reliably 
and obtained an estimate of Vb = 247 ± 19 km s _1 for 
i?o = 7.8 ± 0.4 kpc. Based on a sample of 18 masers, 
Bovy, Hogg & Rix (2009) found V = 244 ± 13 km 



at Rq 



8.2 kpc. Using 44 masers, Bajkova & 



Bobylev (2012) found: (U Q , V Gl W ) = (7.6, 17.8, 8.3) ± 
(1.5, 1.4, 1.2) km s" 1 , Q = -28.8 ± 0.8 km s" 1 kpc" 1 , 
n' = +4.18±0.15kms- 1 kpc- 2 , fi(f = -0.87±0.06km 
s-ikpc- 3 , V Q = \R tt Q \ = 230±14kms- 1 . 

It is important that the rotation-curve parameters found 
are in good agreement with the results of analyzing young 
Galactic disk objects rotating most rapidly around the cen- 
ter: OB associations with SIq = ~ 31 ± 1 km s _1 kpc -1 
(Mel'nik, Dambis, & Rastorguev 2001; Mel'nik & Dambis 
2009), blue supergiants with Q = -29.6 ± 1.6 km s" 1 
kpc- 1 and il' = 4.76 ± 0.32 km s" 1 kpc" 2 (Zabolot- 
skikh, Rastorguev & Dambis 2002), or OB3 stars with 
tt = -31.5 ± 0.9 km s" 1 kpc" 1 , Q' = +4.49 ± 0.12 km 
s" 1 kpc" 2 and fig = -1.05 + 0.38 km s" 1 kpc~ 3 (Bobylev 
& Bajkova 2011). 

The galactocentric radial, V Rrl , and tangential, V$ n (n = 
1, . . . , 44, ) velocities of the masers were determined from 
the relations 



V$ n = U n sin 6 n + (V Q + V n ) cos 9 n , 



(26) 



where U n , V n are the heliocentric space velocities adjusted 
for solar peculiar motion. 

The residual tangential velocities AV$ n are obtained 
from the tangential velocities d26l i minus the smooth rota- 
tion curve that is defined by the Galactic rotation parameters 
fio, ^o; an d ^o f° un d. The radial velocities (l27l i depend 
only on one Galactic parameter Slo and do not depend on 
the rotation curve. As our experience showed, the data are 
so far insufficient to reliably extract the density wave from 
the tangential residual velocities of the masers. Therefore, 
here we determine the spiral density wave parameters only 
from galactocentric radial velocities. 

Figure|2]shows the input radial velocity series Vn n , n = 
1, ...,58, Figure [3] shows transformed velocity series 
V , , n = 1, . . . , 58 and main extracted harmonic. The peri- 
odogram obtained using a modified Fourier transform-based 
method is given in fig. 0] An analysis of this periodogram 
allowed us to estimate the following spiral wave parameters: 
the amplitude of the radial perturbations fji = 7.5 + 1.5 km 
c _1 ; the wavelength (interarm distance in the galactocentric 
direction) A = 2.4 ± 0.4 kpc, and the phase of the Sun in 
the spiral density wave x& = —160 ± 15°. The pitch angle 
of the spiral wave estimated from equation (0 for m = 2 is 
—5.5 ± 1°. The significance level of the peak is p — 0.99. 
Significance level was estimated using the simplest method 
based on Schuster theorem (Vityazev, 2001). The error bars 
are based on a Monte-Carlo simulation of 1000 random re- 
alizations of input data assuming that measurement errors 
obey the normal distribution. Note that the value for /# 
found by fitting the harmonic with A = 2.4 kpc to data is in 
good agreement with relation (I2TI 1. 

Power spectrum, obtained using the GMEM, is shown 
in fig.|5] As we can see, this nonlinear reconstruction algo- 
rithm allowed us to get rid of the side lobes near the main 
peak almost completely and, thus, to increase considerably 
the significance of the extracted periodicity (p = 1.0) with 
A = 2.4 kpc. 

For comparison, in our previous study (Bajkova & 
Bobylev 2012) we have obtained from data on 44 masers 
the following spiral density wave parameters: amplitude 
f R = 7.7 ± 1.7 km s _1 , wavelength A = 2.2*0.4 kpc, 
pitch angle i = —5 ± 0.9°, and the phase of the Sun 
X = -147±17°. 

The parameters of Galactic spiral density wave obtained 
are in good agreement with those found by different au- 
thors by applying different methods to different Galactic 
tracer objects (Mel'nik et al. 2001; Zabolotskikhet al. 2002; 
Bobylev & Bajkova 201 1 and many others). 

6 Conclusions 

We used both Fourier transform-based spectral analysis 
technique and the generalized maximum entropy recon- 
struction method to extract a periodic signal from the 
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Fig. 3 Transformed galactocentric radial velocities vs R 
and extracted main harmonic fitted to the data (solid bold 
line). 




Fig. 4 Periodogram (spectrum power) of masers galacto- 
centric radial velocities. 



galactocentric radial velocities of 58 masers with currently 
known high-precision trigonometric parallaxes, proper mo- 
tions and line-of-site velocities. In accordance with Lin & 
Shu (1964) theory, the extracted periodic signal is associ- 
ated with the Galactic spiral density wave. Masers span a 
wide range of galactocentric distances 3 <R< 14 kpc and 
show a large scatter of position angles 6 in the Galactic XY 
plane, making it necessary an exact accounting for both the 
logarithmic nature of the argument and position angles. As a 
result, we re-determined the main parameters of the Galac- 
tic spiral density wave as follows: amplitude of radial ve- 
locities perturbations /# = 7.5 ± 1.5 km s _1 , wavelength 




Fig. 5 The spectrum power reconstructed by the GMEM. 

A = 2.4 ± 0.4 kpc, pitch angle -5.5 ± 1° and phase of the 
Sun in the density wave x© = —160 ± 15°. 
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